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Conductance of a double quantum dot with correlation-induced 
wave function renormalization 
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Abstract 

The zero-temperature conductance of diatomic molecule, modelled as a correlated double quantum dot attached to noninteract- 
ing leads is investigated. We utilize the Rejec-Ramsak formulas, relating the linear-response conductance to the ground-state 
energy dependence on magnetic flux within the framework of EDABI method, which combines exact diagonalization with ab 
initio calculations. The single-particle basis renormalization leads to a strong particle-hole asymmetry, of the conductance 
spectrum, absent in a standard parametrized model study. We also show, that the coupling to leads V ^ 0.5t (t is the hopping 
integral) may provide the possibility for interatomic distance manipulation due to the molecule instability. 
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Recent progress in nanotechnology made possible to fab- 
ricate small systems attached to leads, and to measure the 
conductance of even a single hydrogen molecule [1]. In such 
systems strong electron correlations can play a decisive role 
[2]. The theoretical understanding of the experiments is far 
from complete, partly because the nanoscale leads geometry 
often introduces a degree of irreproducibility, and partly due 
to many body effects on the intrinsic transport properties 
of such single-molecule devices. In particular, the electron- 
phonon coupling was recently shown shrink the conductance 
peak for a single quantum dot in a Kondo regime [3], since 
the on-site Coulomb repulsion U is effectively renormalized 
in a system ground state. Here we describe a qualitatively 
similar effect in a double quantum dot, produced solely by 
electron-electron interaction only through a completely dif- 
ferent physical mechanism; the correlation-induced wave- 
function renormalization [4]. 

We start with the Hamiltonian which is a generalization 
of the Anderson impurity model [5], and can be written as 



H = Hl + Vl+Hc + Vr + Hr, 



(1) 



where He models the central region, Hi^^r) describes the 
left (right) lead, and V^r) is the coupling between the lead 
and the central region. Both Hi^^r-^ and V^^r) terms have 
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a tight-binding form, with the hopping t and the tunnel- 
ing amplitude V , as depicted schematically in Fig. 1. The 
central-region Hamiltonian He describes a double quantum 
dot with electron-electron interaction, and has a form 

He = {(-a - eVb) X! "' ~ *' X] ('^i<tC2<t + h.c.) 
j=i,2 <T=Ta 

+ U Y^ UiTmi + Knin2 + [Zef/R, (2) 

J=l,2 

where ea is atomic energy, Va is an external gate voltage, 
t' is the internal hopping integral, U and K represents the 
intra- and inter-site Coulomb interactions, respectively, and 
the last term describes the Coulomb repulsion of the two 
ions at the distance R. Here we put Z = 1 and calculate all 
the parameters ea, t' , U, and K as the Slater integrals [6] for 
Is-like hydrogenic orbitals v&is (r) 
where a~ 



\/a^/Tv ex.p{—a\r\), 
^^ is the orbital size (cf. Fig. 1). The parameter 



a is optimized to get a minimal ground state energy for 




Fig. 1. Diatomic molecule as a double quantum dot. 
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Fig. 2. Zero-temperature conductance of the system in Fig. 1 as a function 
of tlie gate voltage Vg and interatomic distance R. Is orbital size a~ is 
optimized variationally. The lead paremeters are t — IRy and V — 0.5i. 

whole the system described by the Hamiltonian (1). Thus, 
following the idea of EDABI method [7], we reduce the 
number of physical parameters of the problem to just a three: 
the interatomic distance R, the gate voltage Va, and the 
lead-molecule coupling V (we put the lead hopping t = 
1 Ry = 13.6 eV to work in the wide-bandwidth limit). 

We discuss now the molecule conductivity calculated 
from the Rejec-Ramsak two-point formula [8]: 



G = Go sin^ TT[E{n) - E{0)]/2A, 



(3) 



where Go ~ 2e^/h is the conductance quantum, A = 
1/N piep) is the average level spacing at Fermi energy, de- 
termined by the density of states in an infinite lead pi^tp), 
E(n) and E{0) are the ground-state energies of the system 
with periodic and antiperiodic boundary conditions, respec- 
tively. E{<j}) is calculated for (p = 0,n within the Rejec- 
Ramsak variational method [8], complemented by the in- 
verse orbital size a optimization, as mentioned above. We 
use typically A'^ = 10^ -^ 10"^ sites to reach the convergence. 

In Fig. 2 we show the conductivity for V = 0.5i, and 
different values of the interatomic distance R. The conduc- 
tance spectrum evolves from the situation well-separated 
peaks corresponding to the independent filling of bonding 
and antibonding molecular orbitals (R < 2ao, where ao is 
the Bohr radius), to the single peak in the intermediate range 
(R ~ 3ao), which decays when t' <^V for large R. 

Probably, the most interesting feature of the spectra de- 
picted in Fig. 2 is their strong asymmetry for small R. 
Namely, the Iow-Vg conductance peak, corresponding to 
the system filling (ni+n2) ~ 1 (one hole) is significantly 
wider than the high-Vc; peak for (ni+n2) ~ 3 (one ex- 
tra electron). Such a particle-hole symmetry breaking is 
the novel feature, observed when including the correlation- 
induced basis optimization, and absent in the parametrized- 
model approach [8,9]. It is also a new feature of a nanosys- 
tem, do not observed in mesoscopic double quantum dots 
[10], where the particle-hole symmetry is perfect. 




GATE VOLTAGE, V'o (Ry/e) 



■2-10 1 

GATE VOLTAGE, Va (Ry/e) 



Fig. 3. The optimized inverse inverse odbital size (a) and average central 
region occupancy (fc) for the diatomic molecule attached to noninteracting 
leads characterized by the parameters t — IRy and V — 0.5t. 

The relation between the observed asymmetry, basis 
renormalization, and electron correlations can be clarified 
as follows. First, one can observe that the optimal values of 
the variational parameter a, provided in Fig. 3a, decrease 
dramatically for high Vg, corresponding to the overdoped 
situation {ni+n2) > 2 (cf. Fig. 3b). This is because the 
system minimize energy of double occupancies, which is of 
the order U ^ a [6]. Then, we focus on an small R limit, 
in which the well separation of molecular orbitals allows 
one to approximate the expression for a single impurity at 
zero temperature [11] 



G = Gosm{7v{nk)/2), 



(4) 



where (uk) is the bonding (k = 0) or antibonding (fc — tt) 
orbital occupancy. Expanding around the maximum we have 



G{Vg) - G{V^-') 



^Go(xcV2)'(Vg 



Vq' ), where 

Vq" (k — 0, tt) is the low/high voltage peak position and 
Xc is the charge susceptibility, which may be approximated 
as Xc ~ d{ni+n2}/dVG, since the orbitals are filled sepa- 
rately. Values of the later derivative read from Fig. 3b around 
the low and high voltage peak positions ((ni+n2) ~ 1 and 
« 3, respectively) provides a clear asymmetry. Moreover, 
the expansion of G(Vg') allows one to roughly estimate the 
peak width as AVg ~ Xc^ ^ {U + K)/2 ~ a, what 
provides an indirect correspondence between the spectrum 
asymmetry and basis renormalization. 

The practical possibility of the conductance measurement 
involving atoms manipulation (changing of R) were also 
explored in terms of the system stability. Namely, the bind- 
ing energies AE = Er^.^ — Eoo (where Rmin is the bond 
length), listed in Table 1, shows that the system become 
metastable around V — 0.5t (AE > indicates a local en- 
Table 1 

The binding energy AE and the bond length i?min for different coupling 
to the lead V . The conductance at the energy minimum is also provided. 



V/t 


AE ii„ 


in/ao 


GlZ 


=°/Go 


n.n 


-0.296 


1.43 




,0-S 


0.1 


-0.293 


1.44 




0.2 


-0.282 


1.46 




0.004 


0.3 


-0.180 


1.52 




0.023 


0.4 


-0.040 


1.59 




0.093 


0.5 


-1-0.074 


1.93 




0.422 



ergy minimum). Tiierefore, tiie binding of the atoms to the 
lead become stronger than interaction between the atoms, 
that may allow individual atoms manipulation. 

In summary, we have supplemented the double quantum 
dot conductance calculation with the single-particle basis 
optimization. Apart from simplifying the analysis, such an 
approach leads to substantially new physical effects, since 
it breaks the particle-hole symmetry of the Hamiltonian. 
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